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ABSTSACT 


The Alpine fault of New Zealand is a major continental transform fault 
which has been uplifted on its southeast side 4 to 11 km within the last 5 
m.y. This uplift has exposed the Haast schists, which have been metamorphosed 
from the adjacent Torlesse graywackes. The Haast schists increase in 
metamorphic grade from prehnite-pumpellyite facies 9—12 km from the fault 
through the chlorite and biotite zones of the greenschist facies to the 
gamet-oligoclase zone of amphibolite facies within 4 km of the fault. These 
metamorphic zone boundaries are subparallel to the fault for 350 km along the 
strike. The K-Ar and Rb-Sr ages of the schists increase with distance from 
the fault; from 4 m.y. within 3 km of the fault to approximately 110 m.y, 20 
km from the fault. Field relations show that the source of heat that produced 
the argon depletion aureole was the fault itself. Adopting a friction model 
for the fault, and applying the known geological history of motion of the 
Alpine fault, we show that the metamorphism resulted from frictional heating 
during the 360 km right lateral slip on the Alpine fault during the Mesozoic 
(Rangitata) period of fault motion (140-120 m.y, ago) , Fault motion began 
again in the Plio-Pleistocene Kaikoura orogeny, with 120 km of further right 
lateral slip and 4-11 km of uplift. Frictional heating during this episode 
produced the argon depletion. Quantitative models for both the metamorphism 
and argon depletion require that the frictional shear stress acting on the 
fault during both episodes of fault motion was at least 1-1,5 kbar. 



INTRODUCTION 


The state of stress within the lithosphere is fundamental to tectonics, 
yet it has remained enigmatic. Whereas the directions of the principal 
stresses at or near the earth's surface can be estimated with a variety of 
geological and geophysical techniques, the magnitude of deviatoric stresses 
that may exist within the lithosphere is as yet uncertain even to their 
order of magnitude. All that is firmly known at present is that shear 
stresses in tectonically active areas must exceed the stress drops of earth- 
quakes, which are mainly confined to the range 10 to 100 bars [Ak'tj 1972; 
Thatohev and Hanks ^ 1973; Kanamovi and And&TSon^ 1975; Hanks ^ 1977]. 

It is particularly important to know the value of the shear stress re- 
quired to produce slip on crustal faults, since the frictional stresses at 
plate boundaries may well provide the primary forces for resisting the mo- 
tion of the lithospheric plates. If these frictional stresses are low, of 
the order of earthquake stress drops, then their resulting resistance to 
plate motion is comparable to several likely driving forces [e.g., Forsyth 
and Uyeda, 1975]. If, on the other hand, they are of the order of kilobars, 
as suggested by laboratory measurements of rock friction [e,g., Sohotz^ 1977; 
ByerteSj 1978] , then friction at plate boundaries becomes a major force re- 
sisting plate motion which may well require that basal shear stresses be an 
important driving force [Hanks^ 1977; Davies ^ 1978]. 

If a fault has slipped a distance u, acting under an average shear 
stress T, then the work per unit fault area, tu, must be dissipated according 
to 

Tu = E + aA + Q 
s 
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where E is the elastic energy radiated by earthquakes, aA is the energy 
s 

required to produce new surfaces of area A and specific surface energy a, 
and Q is heat. Various arguments can be advanced to indicate that the 
first two terms on the right hand side of this energy balance can be ne- 
glected with respect to the third, but we can at least say, conservatively, 
that 

TU < Q . 

If we can measure the frictional heat generated by fault slip, Q, and estimate i 

the fault slip u, we can then obtain a lower bound on x. 

In this study we apply this approach to the Alpine fault of New Zealand, 
one of the major continental transcurrent faults. Within the last 5 m.y,, 

4 to 11 km of uplift on the southeast side of this fault has exposed a sequence 
of rocks (Haast schists) which have a metamorphic history that can be related 
to two episodes of heating due to frictional sliding on the fault. We deve- 
lop the evidence for this interpretation in more detail in the following sec- 
tion, but briefly, our reconstruction is as follows. Between 140 and 120 m.y. 
ago (Rangitata orogeny) some 360 km of right-lateral slip on the Alpine fault 
generated sufficient heat to metamorphose a thick sequence of eugeosynclinal 
grayi<7ackes adjacent to the fault. These metamorphic rocks (the Kaast schists) 
were at that time deeply burled and formed part of a broad ductile shear zone. 

The highest metamorphic grade of the Haast schists, gamet-oligoclase zone 
of amphibolite facies, occurs adjacent to the fault, and the metamorphic grade 
decreases with distance from the fault, through the biotite and chlorite zones 
of the greenschist facies to prehnite-pumpellyite and pumpellyite-actinolite 
subschists that occur 9-12 km from the fault. The fault was reactivated about 
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5 m.y, ago (Kaikoura orogeny) during which a further 120 km of right-lateral 
slip occurred, accompanied by 4 to 11 km of uplift that exposed the Haast 
schists in what are now the Southern Alps. During the Kaikoura period of 
fault motion, frictional heating produced only a slight retrograde metamor- 
phism in the Haast schists but was sufficient to almost completely degas 
these rocks of their radiogenic argon and deplete their micas of strontium 
and thereby reset their K-Ar and Rb-Sr ages to very young values near the 
fault. 

In order to account for both the distribution of metamorphic isograds 
within the Haast schists and the distribution of their K-Ar ages, we have 
quanti(^^tively modeled the frictional thermal events produced by both periods 
of £au>lt motion. We fdLnd that within the permissible range of all variable 
parameters for both periods of fault motion, a minimum average value of 1 to 
1.5 kb for the frictional stress on the Alpine fault, 

GEOLOGICAL HISTORY OP THE HAAST SCHISTS 

The South Island of New Zealand was a site of active subduction during 
the Mesozoic, as indicated by the arc terrane of the New Zealand geos^mcline 
[Landis and Bishop 3 1972]. Tlte arc terrane is preserved in the southern part 
of the South Island (Figure 1) . The Fiordland igneous complex in the west 
consists of a high T/P metamorphic terrane intruded by numerous calc-alkaline 
plutons, and is interpreted as the relic of the Mesozoic arc. To the east of 
this complex is a narrow zone of blueschists and ultramafics that abut a tec- 
tonic suture (Livingstone-Macpherson fault zone) that separates these rocks 
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from the voluminous eugeosynclinal sediments of the Torlessc group. The 
latter have been interpreted as trench sediments [Landis and Bishop^ 1972]. 

At the culmination of this Rangitata orogeny the Otago schists, a 
broad NW-SE striking band of schists derived from the Torlesse group, were 
deformed and uplifted. The Otago schists outcrop in a broad NW-SE striking 
antiform. Their highest metamorphic grade, the biotite zone of the green- 
schist facies is found in the central part of their outcrop area, and their 
metamorphic grade decreases both to the NE and SW from this central zone. 

The Otago schists have K-Ar ages of 110-140 m.y. [Harpev and Landis , 1967; 
Sheppard et 1975] and they probably formed during a compressional event 

which produced shortening normal to the arc near the end of the Sangitata 
orogeny . 

The arc terrane was subsequently offset 480 km by right-lateral slip 
on the Alpine fault iVIellman, 1955; Suggate^y 1963]. We thus see exposed in 
the northwest comer of the South Island a repetition of the entire arc se- 
quence (Figure 1) . There the Marlborough schists lie in the identical struc- 
tural position as the Otago schists, and we interpret them as being their 
offset equivalents. We also notice that a swarm of lamprophyre dikes, dated 
at 120 m.y., has been offset 120 km by the Alpine fault [l^lellman and Cooper, 
1971]. Therefore 360 km of fault motion must have occurred in the Mesozoic, 
between 140 to 120 m.y. B.P. Since the motion of the Alpine fault is pre- 
sently oblique slip, with a component of thrust faulting superimposed on 
right-lateral strike-slip movement [CZark and WeZZman, 1959; Suggate, 1963; 
SohoZz et aZ. , 1973], the remaining 120 km of fault slip is assumed to have 
occurred in the last 5 m.y.. that is, since the initiation of uplift of the 
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Southern Alps in the Kalkoura orogeny [Suggate^ 1963]. The present rate 
of motion on the fault is 5 cra/yr right lateral from sea floor spreading 
data [Christo ffet^ 1971; Hayes and Tatwaniy 1972] to 1 cm/yr from geo- 
logical observations [CXavk and Wetlmany 1959]. Vertical motions may 
be several mm/yr [Suggatey 1963, 1968], 

The Haast (or Alpine) schists are petrographically similar to the 
Otago and Jlarlborough schists in that they are also metamorphic equivalents 
of the Torlesse group sediments. These pelitic schists, however, have a 
distinctly different tectonic and metamorphic history. The Haast schists 
outcrop between the Alpine fault and the main divide of the Southern Alps. 

In the southern part of their outcrop area, where they merge with the Otago 
schists, the folds tighten up and the fold axes are rotated (dextrally) into 
near parallelism with the Alpine fault as they approach the fault [Grindteyy 
1963]. The intensity of deformation of these rocks increases towards the 
fault, and within a few km of the fault have a schistocity and lineation sub- 
parallel to the fault. The style of deformation of the Haast schists is in 
accord with their having been fojrmed under conditions close to pure shear in 
the same sense as that which produced the motion on the Alpine fault, and we 
interpret them as the exposure, at an intermediate structural level, of a 
ductile shear zone associated with the early (Rangitata) period of motion on 
the Alpine fault. Furthermore, unlike the Otago and Marlborough schists, 
the Haast schists were not uplifted until the Kaikoura orogeny. The meta- 
morphic grade of the Haast schists also increases toviards the Alpine fault. 

The metamorphic zone boundaries are subparallel to the fault and their distance 
from the fault varies with the amount of uplift, which is generally greatest 
tov/ards the south, judging from the present elevation of the Alps. In the 
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central Alps, from Arthur’s Pass to Haast Pass, this relationship is the 
most regular (Figure 2) . There garnet-oligoclase zone schists of the lowest 
amphibolite facies occur within a few km of the fault. The schists de- 
crease in grade through the biotite and chlorite zones of the greenschist 
facies to prehnite-pumpeilyite facies subschists which occur 9-12 km from 
the fault. Very close to the fault, isolated pegmatite dikes occur as con- 
cordant lenses in the garnet schists and are interpreted as being due to 
partial melting [Figure 2; Grindleyt 1963]. The K-Ar and Rb-Sr ages of 
the Haast schists also show a regular relationship with the fault. In Fig- 
ure 3, South Island K-Ar ages from Sheppard et al, [1975] are plotted as 
a function of distance from the Alpine fault. The oldest ages for the South 
Island are in the range 140-150 m,y., which probably dates the onset of Ran- 
gitata metamorphism; a second group of ages is in the range 100-120 m,y. 
which probably represents late Rangitata uplift and cooling [Sheppard et at.^ 
1975]. In contrast with the regional ages are a group of very young ages 
from the Haast schists just east of the Alpine fault. The K-Ar whole rock 
and biotite ages from the Haast schists between Arthur’s Pass and Haast Pass 
are shown in Figure 4. These data indicate a strong argon depletion within 
10-15 km of the fault, with a form very similar to that observed in contact 
metamorphic aureoles adjacent to intrusive magmatic bodies [Eartf 1964] , ex- 
cept that the lateral extent of the depletion is much greater. 

These distinctive patterns of metamorphism and argon depletion adjacent 
to the Alpine fault have long been recognized [MasoUj 1961, 1962; Eurley et 
al ,3 1962; Grlndley, 1963; Aronson^ 1965; Egrper and Landis ^ 1967; Sheppard 
et at.:, 1975], They were origina.lly interpreted by mson [1961, 1962], with 
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later elaboration due to Hcrper and Landis [1967], as the result of major 
(20 to 30 km) uplift and rotation/differential uplift of rocks whose meta- 
morphic grade was originally established by their great depth of burial 
in the crust. Since the highest grade (highest temperature) metamorphism 
is closest to the fault, the K-Ar ages can be interpreted as the age at 
which the local temperature dropped beneath the Ar retention temperature 
during uplift, being the most recent for the highest temperature facies 
which require the largest temperature drop and were at the greatest ori- 
ginal depth. 

Sheppard et at. [1975], principally on the basis of the K-Ar ages and 
concentration of ^®Ar, argued that the single stage uplift model described 
above did not entirely explain the observations or otherwise meet separate 
constraints. They proposed a two-stage model in which the metamorphic 
grade of the Haast schists was established between 140 and 120 m.y. ago but 
the argon depletion resulted from rejuvenation of the Alpine fault 4 ± 2 m.y. 
ago. They concluded that the source of the heat necessary to establish the 
observed Ar retention pattern was the fault itself. In support of this two- 
stage model are the following points. 

1. The lamprophyre dikes cut the schists with chilled contacts {Grind- 
teyj 1963], therefore metamorphism of the Haast schists occurred prior to 
120 m.y., whereas the argon depletion occurred after 120 m.y. A two-stage 
process is thus required. 

2. To explain the observed argon depletion spectrum by simple uplift, 
the uplift must have occurred progressively over the past 110 m.y., but no 
schist detritus is found in sediments to the west of the Alpine fault of 
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greater than Pliocene age [Gvi-ndleiji 1963; Suggate^ 1963]. The evidence is 
thus that most of the uplift of the Southern Alps occurred in the last 
5 m.y. 

3. The single stage uplift model requires that the K-Ar age be direct- 
ly related to depth of burial and hence metamorphic grade of the schists. 
Sheppard et at, [1975], however, showed that K-Ar age is solely a function 
of distance from the fault and not of metamorphic grade. This is consistent 
with points 1 and 2 that metamorphism occurred in the Mesozoic and argon de- 
pletion in the Cenozoic. 

/ 1 simple uplift model the temperatures are insufficient to degas 

the in the way observed. Thus Hurtey et at. [1962] calculated the 

argon diffusion parameters for biotite by using the uplift model and obtained 
results some ten orders of magnitude hgiher than that subsequently measured 
directly in the laboratory [Feohtig and Katbttzer^ 1966], Any geothermal 
gradient sufficient to cause argon degassing in the upper part of the section 
would produce rocks of higher metamorphic grade in the lower part of the 
section than those that are actually found [Sheppard et at,^ 1975]. 

5. If the metamorphic facies boundaries were originally horizontal, 
reflecting depth of burial, their present outcrop pattern, which indicates 
steeply dipping facies boundaries, requires a large rotation or differential 
uplift distributed through the schists. But no such rotation on either side 
of the fault is indicated for the lamprophyre dikes, which postdate the 
schists [Grindtey^ 1963; Ivellman and Cooper,^ 1972] and the evidence is nega- 
tive for distributed uplift on shear planes within the schists, since late 
and post-metamorphic folds have not been offset [Gruid-tey, 1963]. 
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A FRICTIONAL HEATING MODEL FOR MESOZOIC METAMORPHISM 
AND LATE CENOZOIC ARGON DEPLETION OF THE HAAST SCHISTS 

Following the argument of Sheppcwd et at. [1975] that the argon deple- 
tion of the Haast schists must have involved frictional heating on the Alpine 
fault during its Plio-Pleistocene (Kaikoura) period of movement, we also note 
that the present metamorphic zone boundaries strongly suggest that the fault 
was a source of heat during the earlier (Rangitata) period of fault motion, 
and that the metamorphism resulted from frictional heating. Otherwise, ac- 
cording to our point 5 above, the present outcrop pattern of these metamor- 
phic rocks is very difficult to explain. 

In Figure 5 we schematically show the history of events in a series of 
vertical sections normal to the Alpine fault. Prior to 140 m,y. ago, we 
show the ancestral position of the Alpine fault in a region of normal geo- 
thermal gradient (Figure 5a). Between 140 and 120 m.y. (Figure 5b), the fault 
slipped 360 km and frictional heating elevated the isotherms in the vicinity 
of the fault, metamorphosing the adjacent rocks. We take 300®C as the preh- 
nite/pumpellyite-greenschists facies boundary, and 450°C as the greenschist- 
amphibolite facies boundary, noting that below several km depth, where the 
metamorphism took place, these transition temperatures are almost independent 
of pressure [Emstj 1975], 

Our models will be constrained so that the 450® isotherm and 30C® iso- 
therm occur at distances L^ and L 2 from the fault at a depth H, where Lj 
and L 2 are the average distances from the fault of the amphibolite-greenschist 
and greenschist-prelinlte/pumpellyite boundaries and H is the Cenozoic uplift 
of the Southern Alps. 
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From 120 to 5 m.y. there was no fault motion, the region cooled to a 
normal geothermal gradient, and the metamorphic facies boundaries were as 
shown in Figure 5c. The Kaikoura movements began 5 m.y. ago, with 120 km 
right-lateral slip and H km of uplift (Figure 5d). We will calculate the 
frictional heating produced by this slip, and using appropriate values of 
the diffusion parameters for argon from biotite, we will require that this 
predict the observed argon depletion profile. This second heating event 
produced some retrograde metamorphism in the schists but did not alter the 
position of the zone boundaries {Grindteyj 1963; She-ppard et dl.^ 1975]. 

Since the metamorphism is independent of the argon depletion, we will further 
constrain these models so that H and t are consistent in modeling both 
processes . 

In the region of the central Southern Alps (Figure 2) the parameters 
Li ® 4-5 km, L 2 = 9-12 km, and H = 4-11 Ion. The lower limit on H is given 
by the present elevation of the Southern Alps, the upper limit is derived 
from our calculations as the greatest depth at which biotite would retain 
radiogenic argon in the interval between the Rangitata and Kaikoura oro- 
genies . 

In the calculations that follow, we allow the parameters H, q^, the 
ambient heat flow, and v, the velocity of fault slip, to vary within the 
allowable limits and calculate the lowest value of shear stress t required 
to produce the raetamorphism and argon depletion. Two models of shear stress 
distribution were tried. In model A (Figure 6a) it was simply assumed that 
a shear stress t is applied to the fault from a depth 0.5 to 15.5 km. In 
the second model, B, it was assumed that the shear stress is governed by a 
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coefficient of friction y so that x = ya where a is the effective normal 

n n 

stress on the fault. In this model then (Figure 6b), x increases linearly 
with depth from 0,5 km to 15.5 km with a mean value, x, at 8 km. Since 
heating from below 15 km did not significantly influence the results we 
assumed that x = 0 at depths greater than 15.5 km. 

In our models the source of heat was assumed to be a sheet of negligible 
thickness at the fault. Although, as described earlier, the Haast schists 
have undergone polyphase deformation that increases in intensity as the 
fault is approached, indicating that some of the transcurrent motion took 
place by penetrative shear, early metamorphic (F2) and post-metamorphic (F3) 
mesoscopic and macroscopic folds are not sheared out except within the 
mylonite zone {Grindley , 1963]. We therefore conclude that the principal 
shearing was accommodated within the narrow (< 1 km thick) zone of mylonites 
at the fault [Reedt 1964] . This agrees with the results of 'luen et al, 

[1978] who showed that the region of intense shear in a viscous shear zone 
is an order of magnitude narrower than the width of the accompanying thermal 
anomaly. Therefore we need not consider shear heating over a broad zone, 

Friational metcanorpkism cotd the Rangitata motion. During the Rangitata mo- 
tion, the fault underwest 360 km of slip in some 20 m.y. In the frictional 
model, if x = ^ ~ [Byevteej 1978], then if is lithostatic, 

we would expect x - 3 kb at 15 km depth, so that in our model B x = 1.5 kb. 
If pore pressures are present , x will be smaller . We accordingly calculated 
heat flow for models in which x ranged from 0.5 to 1.5 kb, and v ranged from 
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2 to 12 cm/yr. For the slower velocities (v < 8 cm/yr) the model reaches 
steady state before the 360 km of motion is complete and an analytical 
solution was used. At higher velocities the temperatures do not reach 
steady state and they have been calculated using a finite difference scheme. 
Both types of calculation are described in detail in the appendix. 

Two representative solutions for the temperature distribution at the 
end of the Rangitata motion are shown in Figure 7. Both began with a back- 
ground geothermal heat flow of 1 HFU, t = 1.5 kb and v = 6 cm/yr. Stress 
model A (Figure 7a) predicts higher temperatures than model B close to the 
fault at shallow depths. The reverse is true at distances greater than a 
few kilometers from the fault. Predicted surface heat flow is also shown. 

The anomaly is concentrated much closer to the fault in model A than in 
model B. An acceptable model iaust show the 450° isotherm 4-5 km from the 
fault and the 300° isotherm 9-12 km from the fault at some depth 11 ^ H 4 km. 
Thus, in the models shown in Figure 7, H = 6-8 km. Variations in the assumed 
value of q^, the initial geothermal heat flow, will not influence these or 
subsequent results other than by compressing or expanding the vertical scale 
in Figure 7, and hence influence values of H, For reasonable values of q^, 
say between 0,8 and 1,2 HFU, this will have a negligible effect on our re- 
sults. The depth H as well as the slip velocity v will therefore be the vari- 
able parameters that will affect the value of t required to produce an accep- 
table temperature distribution. 

In Figure 8 we show t vs. v plots for both model A and model B, Tirhich 
show the values of T and v which yield acceptable temperature profiles at 
the depth H indicated. The greater the value of H, the smaller the value of 
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T required to raise the temperatures to appropriate levels. At low velo- 
cities, the curves follow lines of tv * constant, but at higher velocities 
these curves flatten out and begin to rise at the highest velocities, be- 
cause at high velocities the temperatures do not have time to reach steady 
state even though the same amount of heat, tu, is generated by the fault. 

At velocities higher than about 8 cm/yr, we can find a steady state solu- 
tion (after fault motion has ceased) that produces an acceptable tempera- 
ture distribution, but inevitably these solutions predict transient tempera- 
tures near the fault at depth H higher than the melting point of these rocks 
[750®C, according to vfdttaaej 1976] . Since this is not observed, these high 
velocity models are actually unacceptable. The principal difference between 
model A and model B is that to match the observations , about 1 km greater 
uplift is required for model B for the same values of x and v. 

The shear stress minimum for any model is 500 bars for a profile depth 
of 11 km and is correspondingly higher for shallower depths. Hence we can 
say that the mean effective shear stress on the fault was at least 500 bars. 
In fact we favor models with a velocity less than 'v 8 cm/yr as at higher 
(non-steady state) velocities the width of the amphibolite zone becomes 
large with respect to the greenschist zone, and melting is expected along 
the fault zone at depth H, neither of which is observed. (The rare peg- 
matite dikes indicate that only a very small amount of the lowest melting 
fraction of the schists was melted.) Furthermore, the time available for 
metamorphism becomes very short at the higher velocities. The calctilations 
of the second (Kaikoura) heating episode show that this heating took place 
in too short a time to produce more than a slight retrograde metamorphisra in 


the schists. These considerations lead us to conclude that it is very un- 
likely that T was less than 750 bars. 

Argon devletion and the Kaikoura motion. In discussing the Kaikoura motion 
we use the same models of stress variation and the same range of values, but 
the horizontal velocity v is geologically constrained to be within 2.5 to 4 
cm/yr. Starting with a background heatflow of 1 HFU we solve the heatflow 
equation with one side of the fault moving vertically and being planed off 
by erosion as frictional heating resulting from the horizontal motion takes 
place on the fault. This finite difference calculation is described in the 
appendix. The vertical velocity (W) is chosen so that the depth H of the 
Rangitata metamorphic profile reaches the surface after 5 m,y. or 3 m.y,, 
corresponding to horizontal velocities of 2.5 and 4 cm/yr, respectively. 

Most of our calculations were done using 2.5 cm/yr as the horizontal velo- 
city. A list of the models tried is given in Table 1. 

The temperature history of the profile was followed through the 5 m.y. of 
motion at various distances from the fault. Some representative temperature 
histories are shown in Figure 9. In Figtire 9a, for example, the results are 
for a model with stress distribution A, t = 1.5 kb, moving at 2.5 cm/yr hori- 
zontally and 1.67 mm/yr vertically (H = 8.33 Ion), We plot the temperature 
history of points at various distances from the fault and an original depth of 
8.33 km. The point 1 km from the fault, for example, Is heated from 180°C 
to 275®C after 1 m.y., and subsequently cools. At distances greater than 10 
km from the fault, negligible temperature increases occur. In Figure 9b, the 
results for stress model B are shown with the same values of t and v. For 
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model E, temperatures near the fault do not get as high or last as long as 
in model A, but the opposite is true at greater distances. 

We then use these temperature histories and appropriate values for the 
diffusion coefficient (D 2 o/a^) and activation energy (E) for argon diffusion 
from biotite to calculate the fractional loss of argon (F) since the start 
of motion, (1-F) is directly proportional to K-Ar age for periods as short 
as 100 m.y. [Feohtig and KathttzeT^ 1966] . We calculated argon loss for the 
various models listed in Table 1. For each model we tried a variety of val- 
ues for D 2 o/a^ and E to determine what values were permissible in order to 
fit the data. In Figure 10 we show an example. The calculated curves are 
all for model 1 (Table 1), and we indicate each curve by an upper number, 
log Doo/a^, and a lower nxjmber, E. The family of curves that most closely 
fits the data (solid circles) is within the cross-hatched region, which en- 
closes possible values for the diffusion parameters for this particular model. 
It can be seen that a wide range of D 2 o/a^ and E can be used to fit the rather 
limited data, though D 2 o/a^ and E are not independent variables. We could not, 
however, fit the data with t = 500 bars using any choice of D 2 o/a^ or E be- 
cause for such a low shear stress the temperature history does not vary suf- 
ficiently with distance from the fault. 

With values of t = 1 to 1.5 kb, however, the data could easily be fit. 

In Figure 11 we show the possible values of D 2 o/a^ and E for models assuming 
various values of H, t, and W (Table 1). By possible values we mean that 
they will produce a depletion aureole with almost total argon depletion out to 
3-4 km from the fault and almost no depletion at distances greater than 12 km, 
as the data indicate. Note in Figure 11 that within the range of other para- 
meters, permissible values of D 2 o/a^ and E show a nearly linear relationship 
with one another along the lines of H = constant. This simply reflects that 
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decreasing values of H require greater values for argon diffusivity to 
produce the observed depletion. Direct laboratory measurements of both 
D2 o/^^ and E for argon diffusion from biotite are also plotted (x) in 
Figure 11 [Fveahen and Lippolt,, 1965]. Although in general there are 
problems with applying experimentally determined argon diffusion para- 
meters to geological situations [Musse-btf 1969] at least for phlogopite, 
and TulZis [1977] showed that E for argon diffusion is almost in- 
dependent of pressure and results obtained at high temperatures can be 
extrapolated to temperatures less than 600'’C, If the diffusion parameters 
of Frechen and Lippolt are the most relevant for in situ degassing of argon, 
they lie near the H = 5-6 km curves and hence indicate that the Cenozoic 
uplift of the Southern Alps was 5-6 km, and that a shear stress 7 => 1-1.5 kb 
produced the heat flow required for the argon depletion. This is consistent 
with the earlier calculations since if H = 5-6 km, 1-1,5 kb of average shear 
stress are required to produce the metamorphism, even at the highest slip 
rates (Figure 8) , The fact that Rb-Sr ages show a similar depletion aureole 
near the Alpine fault [Aronson j 1965] reflects the fact that the diffusion 
parameters for Sr diffusion from biotite are similar to those for argon 
[dai’tj 1964; Fedhtig and Kaihitzer^ 1966] and is confirmation that thermal 
degassing is the mechanism responsible for the depletion of '*‘^Ar from these 
rocks, and thermal diffusion of Sr from the mica phases is the mechanism for 
the reduction of Rb-Sr ages of mica concentrates from the Haast schists 
[Aronson^ 1965] . We note that only a very slight retrograde metamorphism 
occurred within 1-2 km of the fault during the Kaikoura orogeny, so that 
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neither bulk recrystalllzatlon nor penetrative deformation occurred in the 
schists during the period of argon depletion [Gpindley 3 1963; Skeppard et at.j 
1975]. These factors give us confidence in the applicability of our simple, 
thermally induced, diffusion calculations. We also take confidence in the 
fact that exactly the same pattern of Rb-Sr and K-Ar depletion are observed 
in contact metamorphic aureoles [Hartj 1964] , 

We conclude, then, that the argon depletion of the schists could not 
have been produced with a shear stress as low as 500 bars. At least 750 bars 
average shear stress is required, and the most likely values, self-consistent 
with experimentally determined values of the diffusivity of argon and with 
the metamorphism, are t = 1-1.5 kb, 

DISCUSSION 

The K-Ar and Rb-Sr age distribution in the Haast schists could only have 
been produced by frictional heating on the Alpine fault during the late 
Cenozoic motion of that fault. An entirely consistent argument shows that 
the metamorphism of these rocks also resulted from frictional heating during 
the Mesozoic period of fault motion. Having established that, we can use 
the argon depletion and metamorphism as independent measures of the frictional 
heating. 

Although it is conceivable that the argon depletion could have been pro- 
duced by mean shear stresses as low as 750 bars acting on the fault, if we 
adopt diffusion parameters for argon anyv;here near consistent with the exper- 
imentally measured values, we must conclude that the shear stresses were at 
least 1-1.5 kbars, and that the original depth of burial of those rocks was 


H - 5-6 km. 
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The metamorphism could have been produced by frictional heating involv- 
ing shear stresses as low as SOO bars, usins the most extreme model permis- 
sible. However, if H = 5-6 km,, then at least 1-1.5 kb of shear stress are 
again required, regardless of the velocity of fault motion. There are also 
several additional arguments that H could not have been as great as 10-11 km. 
If Hwere that great, we should expect that 1) melting would be observed close 
to the fault, which is not observed, 2) the ratio of the width of the amphi- 
bolite zone to that of the greenschist zone would be broader than observed, 

and 3) the volume of sediments derived from the Southern Alps, which presently 
have an elevation of 3-4 km, would be much greater than observed. 

We therefore conclude that it is most unlikely that the mean shear stress 
on the Alpine fault is less than 1-1.5 kb. Even adopting the most extreme 

position, our argtiment Is concerned only as to whether it is 750 bars or 1.5 

kbars. There is no question that it must be at least one order of magnitude 
higher than earthquake stress drops. Furthermore, it is important to point 
out that our estimate is a mirvirrrum one. That is, we have assumed an efficient 
system, all other energy losses, through seismic energy release, creation of 
new surface energy, heat of transformation in metamorphic reactions, and heat 
loss through convection of ground water, have been neglected. For our parti- 
cular problem, we do not believe that the latter is significant, since the 
simple pattern of metamorphic zonation and K-Ar ages with respect to the fault 
demonstrate that conduction is the principal heat transfer process, at least 
at the depths of interest here. 

Our weakest assximption concerning the history of fault motion is that the 
entire 120 km offset of the lamprophyre dikes took place after the onset of 
uplift of the Southern Alps in the Kaikoura orogeny. If the transcurrent 
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component of fault slip began prior to that, we would have to reduce the mean 
slip velocity in our calculations, and consequently increase the value of t 
required to produce the argon depletion. If we assume it took place in a 
shorter period of time than 3 m.y, (Model 7, Table 1) we will again find a 
non-steady state solution which will require a higher value of t. Further- 
more, the actual offset of the dike swarm is somewhat uncertain, owing to the 
difficulty in identifying the boundary of the swarm. Published values range 
from 100 to 160 km; we used the average and most quoted value of 120 km. This 
range of uncertainty, however, will not significantly affect our results. 

Several other observations also suggest, quite independently, the exis- 
tence of high shear stresses on the Alpine fault. Pseudotachylites occur 
commonly within the my Ionite zone in lenses as wide as 3 cm 1976; 

R. H. Sibson, personal communication, 1978]. The presence of these rocks 
probably indicates partial fusion during seismic slip. Although if the par- 
ticle velocities and displacements during a given earthquake are large, shear 
stresses need not be on the order of kilobars to produce Incipient fusion 
[MaKenz’ie and Brune^ 1972; WaZlaae^ 1976], the common occurrence of these 
rocks is generally Indicative of high shear stresses [Sibsorij, 1975]. Further- 
more, preliminary transmission electron microscopic studies of the dislocation 
densities and grain size of the mylonites from the Alpine fault are indica- 
tive of shear stresses on the order of a kilobar (S, IThite, personal communica- 
tion, 1978). 

If our model for and assumptions concerning the frictional metamorphism 
(Rangitata orogeny) and argon-depletion (Kaikoura orogeny) of the Haast 
schists is correct, a substantial asymmetric heat flow anomaly should exist 
over the Alpine fault, the asymmetry being due to the differential uplift of 



the Southern Alps. The absence of a heat flow anomaly localized to the San 
Andreas fault (across a dimension approximately twice the depth of fric- 
tional heating) has led Bvime et at. [1969] and Ladhenbrudk and Sass [1973] 
to conclude that t over the upper 15 to 20 km of the San Andreas fault is 
no more than several hundred bars. This Is a substantially smaller fric- 
tional stress than our minimum estimate for the Alpine fault. On the basis 
of the very similar frictional strengths of crustal rocks (e.g., Byertee^ 
1978], it is difficult for us to accept that the frictional strength of 
crustal fault zones vary by an order of magnitude (or more) , except in those 
cases where (1) nearly lithostatic fluid pressures reduce the effective con- 
fining pressure essentially to zero, possibly the case for the 200 km long, 
aseismically slipping section of the San Andreas fault in central California 
[Inuin and Bonnes, 1975] or (2) the strength of the fault zone is locally 
dominated by low strength, high hydration number clay minerals such as ver- 
miculite and illite, perhaps also the case for the creeping section of the 
San Andreas fault in central California since vermiculite and illite are not 
observed to result in stick-slip. 

There are, however, several additional factors not included in the 
thermal models of Brune et at. [1969] and Laahenbruah and Sass [1973], which 
complicate their interpretation. It is now becoming apparent for much of 
California that the relative motion between the Pacific and Americas plates 
in post-early Miocene times has not been confined to the present day San 
Andreas fault but has been distributed over a series of subparallel faults 
in a zone some 100 km wide. Thus from Point Arguello to Bodega Bay in cen- 
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tral California, extensive slivering of the Salinian block has occurred 
[Johnson and Normarka 1974], with 115 km of right-lateral motion in post- 
early Miocene times for the San Gregorio-Hosgri fault [Gvahcan and Diakin- 
son, 1978] , 50 km of post-Miocene slip on the Rinconada fault [Grahama 
1978] , and suspected slip on the offshore Santa Lucia bank fault [Silver 
et al.a 1978] and Nacimiento fault. Most of the motion on the San Andreas 
fault did not occur until the last 5-6 m.y. [Graham and DiakinsoHa 1978], 
a figure close enough to the thermal time constant such that a well developed 
heat flow anomaly may not yet be recognizable at the surface, or may be 
blurred by thermal anomalies produced by the other active fault strands 
mentioned. A similar, situation exists in southern California where the 
present active strand of the San Andreas fault has probably only been active 
since the opening of the Gulf of California [Atuaterj 1970] and slip is 
distributed on several strands: the Banning-Mission Creek, San Jacinto, 

and Elsinore faults. 

It is well known that a substantial heat flow anomaly exists in the 
Coast Ranges of central California [Laahenbruah and SasSa 1973] but it is 
not confined to the immediate region adjacent to the San Andreas fault. 

The excess heat in this anomaly, which extends 50 km on either side of the 
San Andreas fault, is enough to account for a kilobar or so of frictional 
stress at a displacement rate of 5.5 cm/yr. It is certainly possible that 
the breadth of this anomaly is related, at least in part, to the dispersed 
fault motion described above. Alternatively, as in the case of the unusually 
low heat flow anomalies with respect to those that "should" exist at mid- 
ocean rises for reasonable vertical velocities of mass transport from below 
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[Hanks^ 1971] and the widely accepted explanation of this discrepancy in 
terms of hydrothermal circulation [Anderson and Hobart ^ 1976], ground 
water may also interact with the fault zone heat flow anomaly, at least 
at shallow depth, in such a way as to reduce, broaden, and generally ob- 
scure the localized anomaly calculated on the basis of steady-state, thermal 
conduction models. The coincidence of hot springs with active crustal 
fault zones [Barnes et at.y 1978], while hardly proving this to be the case, 
certainly does not dissuade us from this possibility. Moreover, the deple- 
tion of in metamorphic and plutonic rocks formed within the continental 
crust [TayloVj 19 ] leave little doubt that they have extensively inter- 

acted with ground water. Plainly, however, heat flow measurements in the 
Southern Alps would be of great value in resolving this dichotomy, 

It is obviously of considerable geological importance to ask how common 
is such frictional metamorphism. There is no question that the Alpine fault 
is not an isolated case. Elongated linear zones of my Ionites, many of which 
contain cores of magmatic gneiss, which are concentrically surrounded by 
metamorphic aureoles of just the type studied here have been described in a 
number of other basement terranes, where they are usually referred to as 
deep ductile shear zones [Hioolas et al., 1977]. In fact, the explanation 
of the origin of such features offered by Nicolas et al^ [1977] is virtually 
the same as that employed here, except that they did not have sufficiently 
well defined boundary conditions to solve the problem quantitatively. In 
another example, LeForts [1915], Andrieux et at. [1976] and Hornet and Allegre 
[1976] have argued that frictional heating along the main central thrust of 
the Himalayas produced the inverted metamorphism of that range and produced 
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granites by partial melting. Bivd [1978] calculated the heat generated 
by frictional heating in that case. He used several indirect arguments^ 
to support the assumption that shear stresses along the thrust could 
not exceed about 200 bars and therefore found that Insufficient heat 
would be generated by this mechanism unless the slip velocity on the 
thrust exceeds 30 cm/yr, an unlikely high figure. He therefore suggested 
that some other heat source, such as due to mantle upwelling, must be avail' 
able. We suggest that if the geologic evidence for shear heating of the 
Himalayas becomes compelling that the problem may be resolved if the shear 
stresses on the main central thrust are as high as we find for the Alpine 
fault. This result has in fact been suggested by Graham and England 
[1976] to account for the inverted metamorphism of the Pelona schists be- 
neath the Vincent thrust in California, 

The results we report here are such that it is now important to re- 
examine the structure and metamorphic history of a variety of shear zones 
to determine if shear heating does often play an important role in the dev- 
elopment of such features and if shear stresses roust be generally high 
along such zones. 


The arguments presented by Bird were the paradox first presented by Hubbert 
and Rubey [1959], regarding overthrust faulting, and secondly, the topo- 
graphical argument discussed by Jeffreys [1962] . He differed from Jeffreys 
on the latter, however, since Jeffreys argued that the topographical expres 
sion of the Himalayas requires shear stresses in excess of a kilobar. 
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TABLE 1: Models used for argon depletion calculations. 


Model No. Stress distribution t H Duration of fault 



1 A 1.5 6 5 

2 B 1.5 8 5 

3 A . 1.5 8 5 


u 

H B 1.5 5 5 

5 A 1.5 5 5 

b B 1.5 6 5 

7 B 1.5. 5 3 

a B. 1.0 • 6 5 
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FIGURE CAPTIONS 

Fig. 1 (jeneralized geological map of the South Island of New Zealand. The 
schists of special interest to this study are separated into the broad NW 
trending Otago schists and the, Haast schist groupy outcropping in the 
Southern Alps adjacent to the Alpine fault. Data from N. Z. Geol. Surv. 
(1965) and Landis and Bishop (1962). 

Fig. 2 Geology of the !&ast schist group in the central Southern Alps from 
Arthur's Pass to Haast Pass. 

Fig. 3 All biotite, hornblendey and whole rock K»Ar ages from the South 

Island which reflect Hangitata erogenic events are shown in relation to 
their distance from the Alpine fault (after Sheppard et al. , 1975). 

Fig. 4 Biotite and whole rock K-Ar ages east of the Alpine fault between 
Arthur's Pass and &iast Pass. Data sources: "x", Sheppard et al. . 

(1975); open circles, Hurley et al. . (1962); triangles, Wellman and 
Cooper, (1972); closed circle. Mason (1962); square. Harper and Landis, 
(1967). The 4 m.y. age at 12.5 km fTom the fault is from an isolated 
gneiss body and its significance is suspect (see Sheppard ^ al. , 1973). 

Fig. 5 A sequence of vertical cross sections across the Alpine fault showing 
schematically the development of the present features. A) The initial 
configuration; B) The fault moved 360 km right lateral in the Mesozoic 
and friction caused an elevation of isotherms near the fault. This 
produced the metamorphic pattern shown in C) . In late Cenozoic, D) , 120 
km of right lateral motion and H km of vertical motion produced argon 
depletion of the schists and the present metamorphic outcrop pattern. 
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Fig. 6 The two distributions of shear stress on the fault assumed in the 
calculations. In model A, a shear stress, acts on the fault from 0.5 
to 15.5 km depth. In model B, stress increases linearly with depth from 
0.5 to 15.5 km. 

/ 

Fig. 7 Temperature distributions at the end of the Rangitata period of fault 
motion for stress models A and B. In both cases t = 1.5 kbar and v s $ 
cm/yr. The surface heat flow predicted by the models is also shown. 

Fig. 8 Values of x vs. v at various assumed depths of burial H, required to 
produce the observed pattern of metamorphism: a) results for stress 

model A, b) stress model B. The numbered dots indicate depth, H, for 
var5^ous models which were determined by matching the 300^0 and 450<^C 
isotherms with the present surface metamorphic zone boundaries. The 
determinations of H are subject to uncertainties of about 20!t. The 
crosshatehed bands give the values of H for which the values of t and v 
are appropriate. 

Fig. 9 The temperature history of points at four distances from the fault at 
an original depth of 8.33* km. In these models t = 1.5 kb, v s 2.5 cm/yr 
and w = 1.67 mm/yr. a) Stress model A, b) stress model B. 

Fig. 10 Fractional argon loss as a function of distance frcan the fault for 
model 1 (Table 1) and various values of log (upper number) and E 

(lower number) for diffusion of argon. Dots are observed values, and the 
family of curves within the crosshatehed region are those consistent with 


the data. 
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Fig. 11 Permissible values of the diffusion parameters ^20^^^ ^ 

argon diffusion f^om biotite that fit the observations for the various 
models listed in Table 1. Values of D^q/s^ and E fall along curves of 
constant original depth of burial, H. Direct laboratory values for these 
constants are shown as X (Frechen and Llppolt, 1965). 

Fig. A1 The geometry of the fault zone used for the heat flow calculations. 


APPENDIX 


Al. Finite differenae scheme 

We solve the heet flow equation in a rectangular region whose RHS 
is moving in the z direction with velocity W (W is negative if the RHS 
is uplifted relative to the LHS) . Heat is generated along the boundary 
X * 0, 0 < z < D (Fig. Al) . The equations to be solved are: 

LHS: 

RHS: 

where T is temperature, < is thermal conductivity , p is density and c 
is specific heat. * 

Ve use a simple central -difference finite difference scheme with the 
following boundary conditions: 

z » 0 constant temperature (20*C) 

z « A constant temperature (to give required geothermal 

gradient far from fault) 

X ■ 0 prescribed heat generation q(z) * x(z) * v 


3T 


pc 




|I + w|I.S_72t 

3t 3z pc 


X - ±L 


zero horizontal heat flux 



Typical values for the parameters are: 

L * 40 km A « 30 km grid size « 1 km 

K « 0.005 cal cm”^ s**^ c • 0.2 cal gm”^ 

Y 

, !i 

p * 3 gra cm”* 

1 ram yr”^je W ^ 3 mm yr“*^ 

2.5 cm yr"^^ v ^ 4 cm yr”^ 

The finite difference scheme has been checked by setting the vertical 
velocity, W, to zero and running the scheme to a steady state solution 
(S 10 MY). Thi] solution is in satisfactory agreement with the analytic 
steady state solution (Appendix A2) , the maximum differences being about 

57.. 


A2. Analytic steady state solutions 

The heat flow equation in the area to one side of the fault with 
no vertical velocity is: 

a 2i< a 2<r 

JJ5 + 0 0 < 2 <.A, X > 0 

with boundary conditions 



T - 0 


z • 0, A 


il . -qCz) 

3x 2 k 


X • 0, h/2 < z < D + h/2 



X * 0, 0 < z < h/2, D + h/2 < z < A 


Here q(z) (» t(z) x v) is the heat generated by the fault, so q/2 is the 

heat which flows to each side of the fault, h is the grid spacing used 

in the finite difference scheme and is introduced here so that the analytic solu'- 

tion corresponds as closely as possible to the numerical one with W. = 0, The 

geothermal gradient can be ignored in this calculation as it can simply 

be added to the analytic steady state solution. 

We now look for a solution in the form 


T - 



-nirx/A 

e 


sin 


rnrz 

A * 


For heat flow model A, r(z) 


a S . ' ! Sy?X 

“ .nVK A 


T, and we find 


For heat flow model B, t(z) 



I < z < D 
otherwise 



and we obtain 


I 


2is^ (sin aijP+h/2) . sin IsUlZl) 

n^nVD * A 

- ( (IHh/2) cos - h/2 cos 

ot^kD * * 


A3. Modelling of ergon diffusion 

A particular steady state solution (Appendix A2) for the end of the 
Rangitata orogeny is taken. Then it is run through the finite difference 
scheme (Appendix Al) for 5 MT with a horizontal velocity of 2.5 cm/yr (or 
for 3 MY at 4 cm/yr) with a vertical velocity chosen to bring the re- 
quired metamorphic cross-section to the surface at the end of the motion. 
These velocities are typically 1 to 2 mn/yr. The temperature history 
of this surface is followed, then the temperature history curves are inter 
polated by cubic splines in order to give a sufficiently detailed repre- 
sentation of the temperature for the nxanerical integration below to con- 
verge. About 40 points were required for each curve. The fractional 
argon loss, F, is calculated as 


NAt 

F - SAF • Z 


{-^1^.— - — B) At 
ti ■ 0 v^/2 /t7 


where At is the time step, N is the number of steps, and 


B « — D 
a^ 



,E 293«Tn 
D - D 20 e " '‘R ‘ 293r 

£ • activation energy 
R « gas constant 
D20 

is diffusion parameter at 20®C 

a2 

a is an effective grain size. 


This approximate fbrmula is derived from one given by Feohtig and Kdlbitzev 
[19^&> p. 70] for diffusion from a spherical grain when the initial argon 
concentration is constant throughout the grain. Their formula gave the 
fractional loss after time t at constant temperature T as 

F - - 4 - - T Bt . F < 0.9. 

To obtain our formula we have differentiated with respect to time to give 
the fractional argon loss LY per time At at constant temperature. Then we 
integrate numerically with variable temperature. This should be sufficiently 
accurate provided that the argon loss, AF, within any timestep is « 0.9. 
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